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Abstract. The Parallel Nanowire Detector (PND) is a photon number resolving (PNR) 
detector which uses spatial multiplexing on a subwavelength scale to provide a single electrical output 
proportional to the photon number. The basic structure of the PND is the parallel connection of 
several NbN superconducting nanowires (»100 nm-wide, few nm-thick), folded in a meander pattern. 
PNDs were fabricated on 3-4 nm thick NbN films grown on MgO (T S =400°C) substrates by reactive 
magnetron sputtering in an Ar/N 2 gas mixture. The device performance was characterized in terms of 
speed and sensitivity. PNDs showed a counting rate of 80 MHz and a pulse duration as low as 660ps 
full width at half maximum (FWHM). Building the histograms of the photoresponse peak, no 
multiplication noise buildup is observable. Electrical and optical equivalent models of the device were 
developed in order to study its working principle, define design guidelines, and develop an algorithm 
to estimate the photon number statistics of an unknown light. In particular, the modeling provides 
novel insight of the physical limit to the detection efficiency and to the reset time of these detectors. 
The PND significantly outperforms existing PNR detectors in terms of simplicity, sensitivity, speed, 
and multiplication noise. 

Keywords: superconducting single-photon detector, thin superconducting films, photon number 
resolving detector, multiplication noise, telecom wavelength, NbN, detection efficiency, maximum 
likelihood, estimation. 
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Equation Chapter (Next) Section 1 

1. Introduction 

Photon number resolving (PNR) detectors are required in the fields of quantum communication, 
quantum information processing and of quantum optics for two class of applications. In one case PNR 
detectors are needed to reconstruct the incoming photon number statistics by ensemble measurements. 
This is the case of the characterization of nonclassical light sources such as single photon [1] or «- 
photon [2] state generators or of the detection of pulse splitting attacks in quantum cryptography, 
where an eavesdropper alters the photon statistics of the pulses [3]. In the second case the PNR 
detectors are needed to perform a single-shot measurement of the photon number. Applications of this 
kind are the linear-optics quantum computing [4], long distance quantum communication (which 
requires quantum repeaters [5]) and conditional-state preparation [6]. 

Among the approaches proposed so far to PNR detection, detectors based on charge-integration or 
field-effect transistors [7-9] are affected by long integration times, leading to bandwidths <1 MHz. 
Transition edge sensors (TES [10]) show extremely high (95%) detection efficiencies but they operate 
at 100 mK and show long response times (several hundreds of nanoseconds in the best case). 
Approaches based on photomultipliers (PMTs) [11] and avalanche diodes (APDs), such as the visible 
light photon counter (VLPC) [2, 12], 2D arrays of APDs [13, 14] and time-multiplexed detectors [15, 
16] are not sensitive or are plagued by high dark count rate and long dead times in the 
telecommunication spectral windows. Arrays of SPDs additionally involve complex read-out schemes 
[14] or separate contacts, amplification and discrimination [17]. We recently demonstrated an 
alternative approach [18, 19], the Parallel Nanowire Detector (PND), which uses spatial multiplexing 
on a subwavelength scale to provide a single electrical output proportional to the photon number. The 
device presented significantly outperforms existing PNR detectors in terms of simplicity, sensitivity, 
speed, and multiplication noise. Here we present the working principle of the device (section 2), a 
review of fabrication and of experimental results (section 3 to 5), an extensive analysis of the device 
operation and corresponding design guidelines (section 6) and the first application of a PND to 
reconstruct an unknown incoming photon number statistics (section 7). 

2. Photon Number Resolution principle 

The structure of PNDs is the parallel connection of N superconducting nanowires (N-PND), 
each of which can be connected in series to a resistor R (N-PND-R, figure lb). The detecting element 
is a 4-6 nm thick, 100 nm wide NbN wire folded in a meander pattern. Each section acts as a 
nanowire superconducting single photon detector (SSPD) [20]. In SPPDs, if a superconducting 
nanowire is biased close to its critical current, the absorption of a photon causes the formation of a 
normal barrier across its cross section, so almost all the bias current is pushed to the external circuit. 
In PNDs, the currents from different sections can sum up on the external load, producing an output 
voltage pulse proportional to the number of photons absorbed. 

The time evolution of the device after photon absorption can be simulated using the equivalent circuit 
of Figure la. Each section is modeled as the series connection of a switch which opens on the hotspot 
resistance R hs for a time t hs , simulating the absorption of a photon, of an inductance L kin , accounting 
for kinetic inductance [2 1 ] and of a resistor R . The device is connected through a bias T to the bias 
voltage source V B and to the input resistance of the preamplifier R out . The n firing sections, in pink, all 
carry the same current I f and the N-n still superconducting sections (unfiring), in green, all carry the 
same current I u . I out is the current flowing through R out . 

Let I B be the bias current flowing through each section when the device is in the steady state. If a 
photon reaches the i th nanowire, it will cause the superconducting-normal transition with a probability 
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r|i=T|(I B /Ic (l) ), where n is the current-dependent detection efficiency and l c w is the critical current of the 
nanowire [20] (the nanowires have different critical currents, being differently constricted [22]). 
Because of the sudden increase in the resistance of firing nanowire, its current (l f ) is then redistributed 
between the other N-l unfiling branches and R out . This argument yields that if n sections fire 
simultaneously (in a time interval much shorter than the current relaxation time), part of their currents 
sum up on the external load. 




Figure 1. (a) Circuit equivalent of a N-PND-R. The n firing sections, in pink, all carry the same current I f and the N-n 
still superconducting sections (unfiring), in green, all carry the same current I„. I out is the current flowing through the 
input resistance R ou , of the preamplifier, (b) Scanning electron microscope (SEM) image of a PND with N=6 and 
series resistors (6-PND-R) fabricated on a 4nm thick NbN film on MgO. The nanowire width is w=100 nm, the 
meander fill factor is f=40%. The detector active area is A d =10xl0 um 2 . The devices are contacted through 70nm 
thick Au-Ti pads, patterned as a 50 SI coplanar transmission line. The active nanowires (in color) of the PND-R are 
connected in series with Au-Pd resistors (in yellow). The floating meanders at the two edges of the PND-R pixel are 
included to correct for the proximity effect. 

The device shows PNR capability if the height of the current pulse through R out for n firing 

-(n) -(1) 

stripes I ou t is n times higher than the pulse for one I ou t , i.e. if the leakage current drained by each of 
the unfiring nanowires SI^Iu-Ib is negligible with respect to I B . The leakage current is also 
undesirable because it lowers the signal available for amplification and temporary increases the 
current flowing through the still superconducting (unfiring) sections, eventually driving them normal. 
Consequently, 51^ limits the maximum bias current allowed for the stable operation of the device and 
then the detection efficiencies of the sections. The leakage current depends on the ratio between the 
impedance of a section Z s and R, ut and it can be reduced by engineering the dimensions of the 
nanowire (thus its kinetic inductance) and of the series resistor (see sec. 6). The design without series 
resistors simplifies the fabrication process, but, as Z s is lower, Slik significantly limits the detection 
efficiency of the device. 

3. Fabrication 

NbN films 3-4 nm thick were grown on MgO <100> substrates (substrate temperature 400°C 
[23]) by reactive magnetron sputtering in an argon-nitrogen gas mixture. Using an optimized 
sputtering technique, our NbN samples exhibited a superconducting transition temperature of T c 
=10.5 K for 40-A-thick films. The superconducting transition width was equal to AT C = 0.3 K. 

Both the designs with and without the integrated series resistors were implemented. A 
scanning-electron microscope (SEM) picture of a 6-PND-R fabricated on MgO is shown in figure lb. 
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The size of detector active area (A d ) ranges from 5x5 pin 2 to 10x10 |om 2 with the number of parallel 
branches (N) varying from 4 to 14. The nanowires are 100 nm wide and the filling factor (f) of the 
meander is 40%. The length of each nanowire ranges from 25 to 1 00 pm. 

The three nanolithography steps needed to fabricate the structure have been carried out by 
using an electron beam lithography (EBL) system equipped with a field emission gun (acceleration 
voltage 100 kV, 20 nm resolution). In the first step e-beam lithography is used to define pads 
(patterned as a 50 Q. coplanar transmission line) and alignment markers on a 450nm-thick polymethyl 
methacrylate (PMMA, a positive tone electronic resist) layer. The sample is then coated with a Ti-Au 
film (60 nm Au on lOnm Ti) deposited by e-gun evaporation, which is selectively removed by lift-off 
from un-patterned areas. In the second step, a 160nm thick hydrogen silsesquioxane (HSQ FOX-14, a 
negative tone electronic resist) mask is defined reproducing the meander pattern. The alignment 
between the different layers is performed using the markers deposited in the first lithography step. All 
the unwanted material, i.e. the material not covered by the HSQ mask and the Ti/Au film, is removed 
by using a fluorine based (CHF 3 +SF 6 +Ar ) reactive ion etching (RIE). Finally, with the third step the 
series resistors (85nm AuPd alloy, 50%-each in weight), aligned with the two previous layers, are 
fabricated by lift off via a PMMA stencil mask. Our process is optimized to obtain both an excellent 
alignment between the different e-beam nanolithography steps (error of the order of 1 00 nm) and a 
nanowire with high width uniformity (less than 10% [24]). 

4. Measurement setup 

Electrical and optical characterizations have been performed in a cryogenic probe station with an 
optical window and in a cryogenic dipstick. 

In the cryogenic probe station (Janis) the devices were tested at a temperature T=5 K. Electrical 
contact was realized by a cooled 50 Q microwave probe attached to a micromanipulator, and 
connected by a coaxial line to the room-temperature circuitry. The light was fed to the PNDs through 
a single-mode optical fiber coupled with a long working distance objective, allowing the illumination 
of a single detector. 

In the cryogenic dipstick the devices were tested at 4.2 K. The light was sent through a single-mode 
optical fiber coupled with a short focal length lens, placed far from the plane of the chip in order to 
ensure uniform illumination. The number of incident photons per device area was estimated with an 
error of 5 %. 

The bias current was supplied through the DC port of a 10MHz-4GHz bandwidth bias-T 
connected to a low noise voltage source in series with a bias resistor. The AC port of the bias-T was 
connected to the room-temperature, low-noise amplifiers. The amplified signal was fed either to a 1 
GHz bandwidth single shot oscilloscope or to a 40 GHz bandwidth sampling oscilloscope for time 
resolved measurements and statistical analysis. The devices were optically tested using a fiber - 
pigtailed, gain-switched laser diode at 1.3 \xm wavelength (lOOps-long pulses, repetition rate 26 MHz) 
and a mode-locked Ti:sapphire laser at 700 nm wavelength (40ps-long pulses, repetition rate 80 
MHz). 

Throughout the paper, the single photon detection efficiency of an N-PND (r| ) or of one of its 
sections (n) are defined with respect to the photon flux incident on the area covered by the device 
(active area A d , typically 10x10 |im 2 ) or by one section (A d /N), respectively. 
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5. Device characterization 



Figure 2a shows a single-shot oscilloscope trace of the photoresponse of a 5-PND under laser 
illumination (A.=700 nm, 80 MHz repetition rate). Pulses with five different amplitudes can be 
observed, corresponding to the transition of one to five sections. The measured 80 MHz counting rate 
represents an improvement of three orders of magnitude over most of the PNR detectors at telecom 
wavelength [7, 14, 25], with the only exception of the SSPD array [17]. 

On similar devices, the single-photon detection efficiency (t| ) at A.=1.3 \xm and the dark-counts rate 
DK were measured as a function of the bias current at T= 2.2 K [18]. The lowest DK value measured 
was 0.15 Hz for r\U 2% , yielding a noise equivalent power [26] NEP=4.2xlO" 18 W/Hz 1/2 . This 
sensitivity outperforms most of the other approaches by one-two orders of magnitude (with the only 
exception of transition-edge sensors [25], which require a much lower operating temperature). 

We investigated the temporal response of a 10x10 p.m 2 4-PND-R probed with light at 1.3 |^m 
wavelength using a 40 GHz sampling oscilloscope (figure 2b). All four possible amplitudes can be 
observed. The pulses show a full width at half maximum (FWHM) as low as 660ps. In a traditional 
10x10 yon 2 SSPD, the pulse width would be of the order of 10 ns FWHM, so the recovery of the 
output current I out through the amplifier input resistance is a factor ~4 2 faster (see section 6.2), which 
agrees with results reported by other groups [27, 28]. As shown in section 6.2, the very attractive N 2 
scaling rule for the output pulse duration unfortunately does not apply to the device recovery time. 
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Figure 2 (a) Single-shot oscilloscope trace during photodetection by a 8.6x8 um 5-PND. The device was tested under 
uniform illumination in a cryogenic dipstick dipped in a liquid He bath at 4.2 K. The light pulses at 700 nm form a 
mode-locked Ti: sapphire laser had a repetition rate of 80 MHz. (b) Photoresponse transients taken with a 40 GHz 
sampling oscilloscope while probing a 10x10 um 2 4-PND-R in the cryogenic probe station under illumination with 1.3 
um, lOOps-long pulses from a laser diode, at a repetition rate of 26MHz. The solid curves are guides to the eyes. 



6. PND Design 

We aim at providing a detailed understanding of the device operation and guidelines for the 
design of PNDs with optimized performance in terms of efficiency, speed and sensitivity. 
The first step is to define the relevant parameter space. The width of the nanowire (w=100 nm) and 
the filling factor (f=50%) of the meander are fixed by technology, the thickness of the 
superconducting film (t=4nm) is the optimum value yielding the maximum device efficiency and the 
active area (A d =10 x 10 |im 2 ) is fixed by the size of the core of single mode fibers to which the device 
must be coupled. We consider single-pass geometries (no optical cavity), but the same guidelines can 
be applied to cavity devices with optimized absorption [29]. The parameters of the PND-R that can be 
used as free design variables are: the number of sections in parallel N, the value of the series resistor 
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R and the value of the inductance of each section L . The number of sections in parallel N can be 
chosen within a discrete set of values (N=2, 3, 4, 6, 7, 10, 17), which satisfy the constraints of w, f, 
size of the pixel and that the number of stripes in each sections is to be odd (we consider the geometry 
of Figure lb). The value of L is the sum of the kinetic inductance of each meander L kin and of a series 
inductance which can be eventually added. L kin is not a design parameter, as it is fixed by w, t, f, A d 
and N. If no series inductors are added (bare devices, L =L kin ), the value of L for each N is listed in 
Table 1. 

Table 1. Inductance (L ) and number of squares (SQ) of each section for all possible values of N. The width of the 
nanowires is w=100 nm, the thickness is t=4 nm. The kinetic inductance per square was estimated (L kin /n=90 pH) 
from the time constant of the exponential decay of the output current (T ou t=Tf=L Wn /R „ t , see sec. 6.2) for a standard 
5x5um 2 SSPD [23]. 



N 


Lo 


SQ 


2 


225 nH 


2500 


3 


153 nH 


1700 


4 


117 nH 


1300 


6 


81 nH 


900 


7 


63 nH 


700 


10 


45 nH 


500 


17 


27 nH 


300 



An additional free parameter, relative to the read-out, is the impedance seen by the device on the RF 
section of the circuit R out , which is 50 O (of the matched transmission line) in the actual measurement 
setup (see section 4), but which can be varied in principle from zero to infinite introducing a cold 
preamplifier stage. 

The target performance specifications are the single-photon detection efficiency (n), the signal to 
noise ratio (SNR) and the maximum repetition rate (speed), which must be optimized under the 
constraints that the operation of the device is stable and that it is possible to detect a certain maximum 
number of photons (n^) dependent on the specific application. 
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Figure 3. (a) Simplified circuit of a N-PND-R, where the two sets of n firing and the N-n unfiring sections have been 
substituted by their Thevenin-equivalents. (b-d) Simulated time evolution of I„ (b), I out (c) and I f (d) for a 6-PND-R as 
n increases from 1 to 6. The parameters of the circuit are: L =L Un =81 nH, R =50 fi, R ou ,=50 ft, R h5 =5.5 kft, and 
t hs =250ps. 

A comprehensive description of PND operation should combine thermal and electrical 
modeling of the nanowires [30]. In this work, a purely electrical model (see section 2 and Figure la) 
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has been used to make a reliable guess on how the performance of the device varies when moving in 
the parameter space. 

In this model, the dependence of L kin on the current flowing through the nanowire was disregarded, 
and it was assumed constant. Furthermore, it has been shown [30] that changing the values of the 
kinetic inductance of an SSPD or of a resistor connected in series to it results in a change of the 
hotspot resistance and of its lifetime, eventually causing the device to latch to the normal state. The 
simplified analysis presented here does not take into account these effects, and considers both R hs and 
ths as constant (Ri ls =5.5 k£l, t hs =250ps), and that device cannot latch. However, the results of this 
approach can still quantitatively predict the behavior of the device in the limit where the fastest time 
constant of the circuit x f (see section 6.2) is much higher than the hotspot lifetime (Xf»t hs ), and give a 
reasonable qualitative understanding of the main trends of variation of the performance of faster 
devices (Xf~t hs ). 

In order to gain a better insight on the circuit dynamics (see sec. 6.2) and to reduce the calculation 
time, the N+l mesh circuit of Figure la can be simplified to the three mesh circuit of Figure 3a 
applying the Thevenin theorem on the n firing sections and on the remaining N-n still 
superconducting (unfiring) sections, separately. 

Figure 3b to d show the simulation results for the time evolution of the currents flowing through R out 
and through the unfiring (I u ) and firing (I f ) sections of a PND with 6 sections and integrated resistors 
(6-PND-R) and for the number of firing sections n ranging from 1 to 6. As n increases, the peak 
values of the output current (I out> figure 3b) and of the current through the unfiring sections (I u , figure 
3c) increase. The firing sections experience a large drop in their current (I f , figure 3d), which is 
roughly independent on n. The observed temporal dynamics will be examined in the following 
sections. 

6.1. Current redistribution and efficiency 

— (n) 

Let 8Iik be the peak value of the leakage current drained by each of the still superconducting 
(unfiring) nanowires when n sections fire simultaneously. The stability requirement translates in the 

condition that for each unfiring section I B + olik < I c (as the leakage current increases with n, 

5Iik'™ x) represents the worst case). This limits the bias current and therefore the single-photon 
detection efficiency (n), which, for a certain nanowire geometry (i.e. w, t fixed), is a monotonically 
increasing function of I B /Ic [20]. For instance, in order to detect a single photon (at A,=1.3 |^m, 
T=1.8K) in a section with an efficiency equal to 80% of the maximum value set by absorption (-32%, 

[27]), 8Iik mx) should be made <33% of I B . Therefore the leakage current strongly affects the 
performance of the device and it is to be minimized, which makes it very important to understand its 

— (n) . . 

dependency from the design parameters: 5Iik (N,L ,R ,R out J . 
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Figure 4. Peak value of the leakage current 81 ik drained by each of the still superconducting (unfiring) nanowires (a) 

-(i) 

and of the output current lout (b) when only one section fires plotted as a function of the number of sections in 
parallel N and of the value of the inductance of each section L . The value of the series resistor R and of the output 
resistor R ou , is 50 SI. The orange line highlights bare devices, the colored bars correspond to devices which respect the 
constraints on the geometry of the structure while the grey bars refer to purely theoretical devices which just show 
the general trend. The leakage current and the output current are expressed in % of the bias current I B because they 
are proportional to it. 



The leakage current for n=l is first investigated and its dependency on n is then presented for 

some particular combinations of design parameters. The dependency of §lL' on N and L at fixed R 
and R out (both equal to 50 O) is shown in figure 4a: an orange line highlights bare devices (L =L kin , 
see Table 1) and the colored bars are relative to devices which respect the constraints on the geometry 
of the structure (L >L kin ), while the grey bars refer to purely theoretical devices which just show the 
general trend. For any N, the current redistribution increases with decreasing L , as the impedance of 

—(i) 

each section decreases. Keeping L constant, Slik decreases with increasing N, as the current to be 
redistributed is fixed and the number of channels draining current increases. For this reason also the 
increase of redistribution with decreasing L becomes weaker for high N. 
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The dependency of 5Iik on R is shown in figure 5 a for some bare devices and R out =50 Q. As 
expected, the redistribution decreases as Ro increases because the impedance of each section increases 

with respect to the output resistance. For the same reason, 5Iik is strongly reduced (to -3% of I B ) 
when Ron, is decreased of one order of magnitude from 50 to 5 Q, keeping R constant (figure 6a). 
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Figure 6. Variation of the peak value of the leakage current per unfiring section (a) and of the output current (b) of 

—(i) -(0 

the set of bare devices for n=l (8Ii k and I ouf , respectively) as the resistance of the output resistor R 0llt decreases of 
one order of magnitude from 50 to 5 fl (in blue and orange, respectively), while R =50 fi. 



— (n) 

The variation of the leakage current with the number of firing stripes n ( 8Iik ) for the set of 

bare devices is presented in figure 7a. The dependency is superlinear (8lL^>n8lL^), as the current to 
be redistributed per firing stripe is always the same (see sec. 6.3), but the number of channels draining 
current decreases. Furthermore, as expected, the curves for different design parameter sets never 
cross, which means that all the design guidelines presented in Figure 5a, Figure 6a, Figure 7a for n=l 
still apply for higher n. 



11 



Physics and application of PNR detectors based on superconducting parallel nanowires 




Figure 7. Variation of the leakage current 81 i k (a) and of the output current I„„i (b) with the number n of firing 
stripes for the set of bare devices. 

In conclusion, the result of this simplified analysis is that, in order to minimize the leakage current 
and thus maximize the efficiency, N, L and R must be made as high as possible and R out as low as 
possible. We note however that R cannot be increased indefinitely to avoid that the nanowire latches 
to the hotspot plateau before I B reaches I c [23]. 

6.2. Transient response and speed 

Before proceeding to the analysis of the SNR and speed performances of the device, it is 
necessary to discuss the characteristic recovery times of the currents in the circuit. 
The transient response of the simplified equivalent electrical circuit of the N-PND (figure 3a) to an 
excitation produced in the firing branch can be easily found analytically. Therefore, the transient 
response of the current through the firing sections If, through the untiring sections I u and through the 
output I out after the nanowires become superconducting again (t>t hs ) can be written as: 

N-n / . \ n I i \ 
I f oc — — exp ( -t I r ) + — exp ( -t I x ; ) 

< I u ocexp(-£/r )-exp(-f/rl (1) 

1 a« I T f ) 

where x s =L /Ro and Tf=L /(R +NR out ) are the "slow" and the "fast" time constant of the circuit, 
respectively. 

This set of equations describes quantitatively the time evolution of the currents after the healing of the 
hotspot in the case if»t hs , and it provides a qualitative understanding of the recovery dynamics of the 
circuit for shorter x f . 

The recovery transients (t>t hs ) of I out , Slik and I f for a 4-PND-R simulated with the circuit of figure 3a 
are shown in figure 8a, b, c, respectively (in blue) for different number of firing sections (n=l to 4). 
As n increases from 1 to 4, the recoveries of I out and 8Ii k change only by a scale factor. On the other 
hand, the transient of If depends on n and becomes faster increasing n, as qualitatively predicted by 
the first of equations (1). Indeed, If consists in the sum of a slow and a fast contribution, whose 
balance is controlled by the number of firing sections n. To prove the quantitative agreement with the 
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analytical model in the limit Tf»t hs , the simulated transients of I out , 8Ii k and I f have been fitted (figure 
8a, b, c, respectively, in orange) using the set of equations (1), and four fitting parameters (x s , if, a 
time offset t and a scaling factor K). The values of x s and if obtained from the three fittings (of I out , of 
Slit and of the whole set of four I f for n=l,..., 4) closely agree with the values calculated from the 
analytical expressions presented above and the parameters of the circuit (x s =2.30 ns ,x f =460 ps). 
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Figure 8. Recovery transients (t>t hs ) of I out (a), 8I| k (b), and I f (c) for a 4-PND-R as n increases from 1 to 4. The 
simulated transients are in blue, the fitted curves are in orange. The parameters of the circuit used for the 
simulations are: L =L kin =117 nH, R =50 ft, R out =50 ft, R hs =5.5 kft, and t hs =250ps. The three sets of curves are fitted 
by equations (1) (multiplied by K, and shifted by t ), where the values of x s and x f are shown in the insets. 



In order to quantify the speed of the device, we take fo=(t r eset)" 1 as the maximum repetition 
frequency, where t reset is the time that I f needs to recover to 95% of the bias current after a detection 
event. 

According to the results presented above, which are in good agreement with experimental data (figure 
2b), I out decays exponentially with the same time constant for any n (x 0U =Xf), which, for a bare N- 
PND, is N 2 times shorter than a normal SSPD of the same surface [27, 28]. This however does not 
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relate with the speed of the device. Indeed, t reset is the time that the current through the firing sections 
If needs to rise back to its steady-state value (IfT b ). In the best case of n=N, If rises with the fast time 
constant x f , but in all other cases the slow contribution becomes more important as n decreases (see 
figure 3d and figure 8c), until, for n=l, I f ~[l-exp(-t/x s )]. The speed performance of the device is then 
limited by the slow time constant (t reset ~3-x s ), which means that an N-PND is only N times faster than 
a normal SSPD of the same surface, being as fast as a normal SSPD whose kinetic inductance is the 
same as one of the N section of the N-PND. 

Figure 9 shows the dependence of fo on L and Ro. For x s <t hs (i.e. fo>4 GHz in our model) the speed of 
the device may be limited by the hotspot temporal dynamics, and so no reliable predictions can be 
made using our simplified model. 




Figure 9. Dependence of f on L and R . No data are presented for f >4 GHz, where no reliable predictions can be 
made using this simplified model. 



6.3. Signal to noise ratio 

The peak value and the duration of the output current pulse are a function of the design 
parameters (see below and section 6.2, respectively). As the output pulse becomes faster, amplifiers 
with larger bandwidth are required and thus electrical noise become more important. In order to assess 
the possibility to discriminate the output pulse from the noise, we define the signal to noise ratio 

(SNR) as the ratio between the maximum of the output current I ou t and the rms value of the noise- 
current at the preamplifier input I n , SNR=I ou t /I n . 

The peak value of the output current when n sections fire simultaneously (see figure 3b, 
relative to a 6-PND-R) can be written as: 

I { :l^{l,-lf)-(N-n)SI^ 

where the starred values refer to the time t=t when the output current peaks. 
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As n=l represents the worst case, in order to evaluate the performance of the device in terms of the 

-(i) -(i) , . 

SNR, the dependency of I ou t from the design parameters is investigated: lout (N,L ,R ,R out J . The 

-(i) 

dependency of I ou t on N and L at fixed Ro and R ou t (both equal to 50 £2) is shown in figure 4b. 

-(i) —(i) 

Inspecting the values of I ou t and of Slik for the same device in figure 4, it becomes clear that they 
add up to a value well above to 1 B , which is due to the fact that the output current and of the leakage 
current peak at two different times, t* and t lk , respectively (figure 3). Furthermore, as tn>t*, the output 
current is not significantly affected by redistribution, because l out is maximum when 8Ii k is still 
beginning to rise. 

The expression for tik can be derived from eq. (1): ti k =L /(N-R out )ln(l+N-R ou t/Ro), which means that 
increasing the device speed (decreasing L or R ), N or R out makes the redistribution faster and then 

lout lower. 

So, for any given N, iii decreases (figure 4b) with decreasing L , both because 8lL' is higher and 

because t [k is lower. Keeping L constant, iLi decreases with increasing N because even though SlL' 
decreases, the redistribution peaks earlier and the number of channels draining current increases. 
-(i) 

The dependency of I ou t on R is shown in figure 5b for some bare devices and R ou t = 50 £1 Even 

though Slik decreases as R increases (figure 5a), the output current decreases due to the 
redistribution speed-up (decrease of t lk ): 8Ii k (1) * increases despite of the decrease of the peak value of 
the leakage current. On the other hand, a decrease in R out makes the redistribution much less effective, 
as tik decreases slower with decreasing R out than with increasing R . Indeed, as shown in figure 6b for 

-(i) 

bare devices, I ou t significantly increases when R out is decreased by one order of magnitude from 50 to 
5 Q, keeping R constant. 

In conclusion, in order to maximize the output current, N, R and R out must be minimized, while L 
must be made as high as possible. 

The rms value of noise-current at the preamplifier input l n can be written as I n =^S n Af , 

where S n is the noise spectral power density of the preamplifier and Af is the bandwidth of the output 
current I out , which is estimated as Af=l/x out , where T ou t = T^Lo/(Ro+NR out ) is the time constant of the 
exponential decay of I ou t (see sec. 6.2). I n is then a function of the parameters of the device and of the 
read-out through S n and x f . As S n varies consistently with the type of preamplifier used, in our analysis 
of I n we take into account only the variation of Af. Assuming S n constant, I n is thus minimized 
minimizing N, R and R out and maximizing L . 

The same optimization criteria apply then naturally to the SNR. The dependence of I ou t/"v/Af on N 
and L for an input resistance of Rout=50£2 is shown in figure 10. 
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The main design guidelines which can be deduced from the analysis of sections 6.1 to 6.3 are 

summarized in Table 2. The type of dependency of 5Iik , fo, lout and Af from the design parameters 
(Lo, R-o, Rout, N) is indicated. 

Table 2. Dependency of 8Iik , fo, lout and Af from the design parameters: increasing with increasing the 
parameter (71), decreasing with increasing the parameter ("ii), independent ( — ). 





Lo 


Ro 


Rout 


N 


Slik 


"^l 


"^l 


71 




fo 


"^l 


71 






lout 


71 


"^1 


■^1 


"^i 


Af 




71 


71 


71 



7. Application to the measurement of photon number statistics 

We wish to determine whether the PND can be used to measure an unknown photon number 
probability distribution. Indeed, the light statistics measured with a PND differ from the original one 
due to non-idealities such as the limited number of sections and limited and non-uniform efficiencies 
(r|i) of the different sections. As a PND can be modeled as a balanced lossless N-port beam splitter, 
every channel terminating with a single photon detector [19], the input-output transformation can be 
formalized as follows. 

Let an N-PND be probed with a light whose photon number probability distribution is 
S=[S(m)]=[s m ], and its output be sampled H times. The result of the observation can be of N+l 
different types (i.e. 0,..., N stripes firing), so an histogram of the H events can be constructed, which 
can be represented by a (N+l)-dimensional vector r=[rj, where r ; is the number of runs in which the 
outcome was of the i th type. The expectation value of the statistics obtained from the histogram is 
E[Q ex =r/H]=Q, where Q=[Q(n)]=[q n ] is the probability distribution of the number of measured 
photons. 

Q(n) is related to the incoming distribution S(m) by the relation: 
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Q(n) = ^P N (n\m)-S(m) (2) 

in 

where P N (n|m) is the probability that n photons are detected when m are sent to the device. Eq. (2) 
may be rewritten in a matrix form as Q=P N S, where p N = j^p N (n|m)J = ^p^ 11 J is the matrix of the 
conditional probabilities. 

It has been shown [31, 32] that an unknown incoming photon number distribution S can be recovered 
if Q and P N are known. 

7.1. Matrix of conditional probabilities 

The matrix of the conditional probabilities of a N-PND depends only on the vector of the N 
single-photon detection efficiencies of the different sections of the device T|=[r|i] through the relations 
presented in [18, 19]. The vector r\ can be then determined from the statistics Q ex measured when 
probing the device with a light of known statistics S. 

For example, using a laser light source with Poissonian photon number probability 
distribution, the probability distribution of the number of measured photons Q (expressed by (2)) was 
fitted to the experimentally measured distribution Q ex using r\ as a free parameter. The resulting r\ and 
matrix of the conditional probabilities are shown in figure 1 1 for a 5-PND at ^.=700 nm. 




Figure 11. Conditional probability matrix for a 8.6x8 urn 2 5-PND (with no integrated series resistors), calculated from 
the vector r\ of the 5 single-photon detection efficiencies (relative to T=4.2 K, ^=700nm) of the different sections of the 
device (inset). 

The P N matrix provides a full description of the detector. Once P N is known, several approaches can 
be used to reconstruct S from the histogram r. In the case no assumptions on the form of S are made, 
the maximum likelihood (ML) method is the most suitable, as it is the most efficient in solving this 
class of problems [33]. 
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7.2. Maximum-Likelihood (ML) method 

Let R=R ,..., Pv N be the random vector of the populations of the (N+l) different bins of the 
histogram after H observations. The joint probability density function L(r|Q) for the occurrence of the 
particular configuration r=r ,. . ., r N of R is called the likelihood function of r and it is given by [34]: 



L(r|Q) = /nnff 



(3) 



1=0 



where Q=[q;] is the probability distribution of the number of measured photons, i.e. the vector of the 
probabilities to have an outcome in the bin i (i=0,. . ., N) in a single trial. 

Considering Q as a function of S through equation (2), we can rewrite the likelihood function of the 
vector r, depending on the parameter S: 



L(r|S) = ^n 



N s 

m m 



V m 



(4) 



i=0 



which is then the probability of the occurrence of the particular histogram r when the incoming light 
has a certain statistics S. 

As r is measured and then it is known, L(r|S) can be regarded as a function of S only, i.e. L(r|S) is the 
probability that a certain vector S is the incoming probability distribution when the histogram r is 
measured. The best estimate of the incoming statistics which produced the histogram r according to 
the ML method is the vector S e which maximizes L(r|S), where r is treated as fixed. So, the 
estimation problem can in the end be reduced to a maximization problem. 

For numerical calculations, it is necessary to limit the maximum number of incoming photons 
to m max (in the following calculations, m max =21). As S is a vector of probabilities, the maximization 
must be carried out under the constraints that the s n are positive and that they add up to one. The 

2 

positivity constraint can be satisfied changing variables: S n -0~ n . Instead of L, we maximize the 
logarithm of L: 



7(E) = ^(Z)) = ln(C) + 5>,ln |>,>, 



\m=0 



(5) 



where £=[o n ] an d C is a constant. 

The condition that the s n add up to one can be taken into account using the Lagrange multipliers 
method: F(S,a)=l(S) -a £o^-l 



\ m=0 



After developing [35] the set of m max +2 gradient equations VF(E,a)=0, we obtain that a=H and 
that the set of m max +l nonlinear equations to be solved respect to I is: 
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y^Pu H 

i—t m „,ax 

/ i Jri,m m 
m=0 



(6) 



for 1=0,. . ., m max . The set of equations (6) can be solved by standard numerical methods. 



7.3. ML reconstruction 

To test the effectiveness of the reconstruction algorithm, a 8.6x8 |^m 2 5-PND was tested with the 
coherent emission from a mode -locked Ti:sapphire laser, whose photon number probability 
distribution is approximately poissonian (S(m)=|a m -exp(-|j.)/m!). Therefore, S could be determined by 
measuring the mean photon number per pulse \i with a powermeter. To determine Q ex , histograms of 
the photoresponse voltage peak were built for varying |i. The signal from the device was sent to the 1 
GHz oscilloscope, which was triggered by the synchronization generated by the laser unit. The 
photoresponse was sampled for a bin time of 5ps, making the effect of dark counts negligible. 
The device was characterized in terms of its conditional probability matrix P 5 ([18, 19], figure 1 1), so 
it was possible to carry out the ML estimation of the different incoming distributions with which the 
device was probed. Because of the bound on the number of incoming photons which is possible to 
represent in our algorithm (m max =21) and as for a coherent state losses simply reduce the mean of the 
distribution, the ML estimation was performed considering |i*=|j./10 and T]*=l Oil (the efficiency of 
each section being lower than 10%). 

Figure 12 shows the experimental probability distribution of the number of measured photons 
Q ex obtained from the histograms measured when the incoming mean photon number is |a=1.5, 2.8, 
4.3 photons/pulse (figure 12 a, b, c respectively, in orange), from which the incoming photon number 
distribution is reconstructed. The ML estimate of the incoming probability distribution S e is plotted in 
figure 12 d, e, f, (light blue), where it is compared to the real incoming probability distribution S 
(green). The estimation is successful only for low photon fluxes (|i=1.5, 2.8 figure 12 d, e) and it fails 
already for u=4.3 (figure 12 f). In figure 12 a, b, c, Q ex (orange) is compared to the ones obtained 
from S and S e through relation (2) (Q, Q e in green and light blue, respectively). 
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Figure 12. (a), (b), (c) Probability distribution of the number of measured photons obtained from experimental data 
Q ex (orange), from S (Q, in green) and S e (Q e5 in light blue) through relation (2) for u.=1.5, 2.8, 4.3 photons/pulse, 
respectively, (d), (e), (f) Real incoming probability distribution S (green) and its ML estimate Se (light blue) for u.=1.5, 
2.8, 4.3 photons/pulse, respectively. The 8.6x8 um2 5-PND was tested under uniform illumination in a cryogenic 
dipstick dipped in a liquid He bath at 4.2 K. The light pulses at 700nm form a mode-locked Tirsapphire laser were 
40ps wide (after the propagation in the optical fiber) and the repetition rate was 80MHz. The average input photon 
number per pulse u was set with a free space variable optical attenuator. 



The main reasons why the reconstruction fails are not the low efficiencies of the sections of 
the PND or the the spread in their values, but rather the limited counting capability (N=5) and a poor 
calibration of the detector, i.e. an imperfect knowledge of its real matrix of conditional probabilities. 
This assessment is supported by the following argument. If we generate Q ex with a Montecarlo 
simulation [19] using the same r\ vector of Figure 11 and a poissonian or thermal incoming photon 
number distributions and then we run the ML reconstruction algorithm (using the same P 5 , which this 
time describes perfectly the detector), S can be estimated up to much higher mean photon numbers. 
However, to alleviate this problem, a self-referencing measurement technique might be used [36]. 
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8. Conclusion 

A new PNR detector, the Parallel Nanowire Detector, has recently proven to significantly 
outperform existing approaches in terms of sensitivity (NEP=4.2xlO" 18 W/Hz 12 ), speed (80 MHz 
count rate) and multiplication noise [18, 19] in the telecommunication wavelength range. 
An electrical equivalent circuit of the device was developed in order to study its operation and to 
perform its design. In particular, we found that the leakage current significantly affects only the PND 
detection efficiency, while it has a marginal effect on its signal to noise ratio. Furthermore, in order to 
gain a better insight on the device dynamics, the (N+l)-mesh equivalent circuit of the N-PND was 
simplified and reduced to a three mesh circuit, so that the analytical expression of its transient 
response could be easily found. With this approach, we could predict a physical limit to the recovery 
time of the PND, which is slower than that previously estimated. Additionally, the figures of merit of 

the device performance in terms of efficiency, speed and sensitivity ( 5Iik , fo, SNR) were defined and 
their dependency on the design parameters (L , Ro, R ou t, N) was analyzed. 

In order to prove the suitability of the PND to reconstruct an unknown light statistics by ensemble 
measurements, we developed a maximum likelihood estimation algorithm. Testing a 5-PND with a 
Poissonian light we found that the reconstruction of the incoming photon number probability 
distribution to be successful only for low photon fluxes, most likely due to the limited counting 
capability (N=5) and the poor calibration (i.e. the imperfect knowledge of the real matrix of 

conditional probabilities) of the detector used, and not to its low detection efficiency (nD3%). 

Additional simulations will be needed to evaluate the performance of our detector for the 
measurement of other, nonclassical photon number distributions. Finally, despite the high sensitivity 
and speed of PNDs, their present performance in terms of detection efficiency (n=2% at A,=1.3|im) 
does not allow their application to single shot measurements, as required for linear-optics quantum 
computing [4], quantum repeaters [5] and conditional-state preparation [6]. Nevertheless, the n of 
SPDs based on the same detection mechanism can be increased to -60% [29], and could potentially 
exceed 90% using optimized optical cavities. 
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